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We present a study of the decay B~^D°(yp^K~ and its charge conjugate, where ^"cp) is recon- 
structed in both a non-CP flavor cigcnstatc and in CP (CP-cvcn and CP-odd) eigenstates, based on 
a sample of 382 million T{4S) BB decays collected with the BABAR. detector at the PEP-11 e+e" 
storage ring. We measure the direct CP asymmetries Acp± and the ratios of the branching frac- 
tions Rcp±: Acp+ = 0.27±0.09(stat) ±0.04(syst), Acp- = -0.09±0.09(stat) ±0.02(syst), Rcp+ = 
1.06 ± O.lO(stat) ± 0.05(syst), Rcp- = 1.03 ± O.lO(stat) ± 0.05(syst). We also express the results in 
terms of the so called Cartesian coordinates x+, x-, and r^: x+ = —0.09 ± 0.05(stat) ± 0.02(syst), 
X- = 0.10 ± 0.05(stat) ± 0.03(syst), r^ = 0.05 ± 0.07(stat) ± 0.03(syst). These results will help to 
better constrain the phase parameter 7 = aig{—VudV*ii/VcdV*f,) of the Cabibbo-Kobayashi-Maskawa 
quark mixing matrix. 
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The angle 7 = a.i:g{-VudV*i,/VcdV*^) is one of the least 
precisely known parameters of the corresponding uni- 
tarity triangle of the Cabibbo-Kobayashi-Maskawa ma- 
trix [1]. There are many proposals on how to measure 7 
involving charged B decays. The B~^D^*^^K^*^~ de- 
cay mode [2], which exploits the interference between 
b — > CMS and b ucs decay amplitudes, is one of the 
most important of these [3, 4]. In this paper we use a 
theoretically clean measurement technique suggested by 
Gronau, London, and Wyler (GLW). It exploits the in- 
terference between B~^D^K~ and B~^D^K~ decay 
amplitudes, where the and mesons decay to the 
same CP eigenstate [3]. We express the results in terms 
of the commonly used ratios Rcp± of charge-averaged 
partial rates and of the partial-rate charge asymmetries 

^CP±, 



Rcp± 
Acp± 



T{D--.D[\^, . K-) + T{B+-^D[\^, , K+) 



[r{B-^D"K 
T{B 



D"cP±K- 



r{B+^D"K+)] /2 
- T{B+ 



D"cP±K+) 



T{B- 



^D%P±K- 



^D'cp±K+) 



(1) 



(2) 



Here, D^p^ 



(D" ± D°)/V2 are the CP cigen- 
states of the neutral D meson system, following the 
notation in Ref. [5]. Neglecting D^—D'^ mixing [6], 
the observables Rcp± and Acp± are related to the an- 
gle 7, the magnitude ratio r of the amplitudes for the 
processes B~^D'^K~ and B~^D^K~ , and the rela- 
tive strong phase 6 of these amplitudes through the re- 



lations Rcp± = 1 



± 2r cos (5 cos 7 and Acp± = 



±2r sin^ sin7/i?c'P± [3]. Theoretical predictions for r 
are on the order of 0.1 [3], in agreement with recent 
results by BaBar (r = 0.091 ± 0.059 [7]) and Belle 
(r = 0.159 ± 0.074 [8]), obtained through the study of 
B-^D^K-, £)0 K+TT-TT^ and D° ^ KItx+tt- de- 
cays. 

This analysis, based on 348 fb"""" of data collected at 
the T{AS) resonance, updates a previous BABMi study 
based on 211 fb~^ of data [9]. Belle recently presented 
a similar measurement of Rcp± and Acp± based on 
251 fb"^ of data [10]. 

The ratios Rcp± are computed under the assumption 
Rcp± = R±/R. which holds neglecting a factor of < 
0.012 as discussed later. The quantities R+, R-, and R 
are defined as: 



(±) 



BiB- 



^Dlcp±)K-)+BiB+^Dj^p^^K+) 
-^Dl.p^^7r-)+B{B+^Dlp^^7r+) 



(3) 

Several systematic uncertainties affect the D^K and D'^n 
final states in the same way and therefore cancel in the 

double ratios Rcp+ and Rcp- , for instance the uncertain- 
ties on charged particle reconstruction efficiencies, and 



the uncertainties on the secondary branching ratios of 
the decays. We express the CP-sensitive observables 
in terms of three independent quantities x+, X-, and r: 

Rcp+{1 T Acp+) - -Rcp-(1 T Acp-) , . 



2 2,2 RCP+ + RCP- - 2 

r = x±+y± = 



(5) 



where x± = r cos((5 ± 7) and y± = r s'm{S ± 7) are the 
so called Cartesian coordinates related to the CP pa- 
rameters that are measured using a Dalitz analysis of 
B-^D^K-, £)0 ^ 7^:077-77+ decays [8, 11]. This choice 
allows the results of the two measurements to be ex- 
pressed in a consistent manner. 

The measurements use a sample of 382 million T(4S') 
decays into BB pairs collected with the BaBar de- 
tector [12] at the PEP-II asymmetric-energy B fac- 
tory. Chargcd-particle tracking is provided by a five- 
layer double-sided silicon vertex tracker and a 40-layer 
drift chamber (DCH). A ring- imaging Cherenkov de- 
tector (DIRC) provides additional particle identification 
(PID). Photons are identified by the electromagnetic 
calorimeter (EMC), which is comprised of 6580 thallium- 
doped Csl crystals. These systems are mounted inside 
a 1.5T solenoidal superconducting magnet. Wc use the 
GEANT [13] software to simulate interactions of particles 
traversing the detector, taking into account the varying 
accelerator and detector conditions. 

We reconstruct B~—^D^h~ decays, where the prompt 
track h~ is either a kaon or a pion. The candi- 
dates arc reconstructed in the CP-even cigcnstatcs tt^tt"*" 
and K^K+ (D^p+), in the CP-odd eigenstates if°7r° 
and K2;Uj {D^p_), and in the (non-CP) flavor eigen- 
state A'^7r+. The lo candidates are reconstructed in the 
7r~7r"'"7r° channel, and candidates in the tt+tt" chan- 
nel. Compared to the previous analysis [9], the current 
study does not include the decay mode Kl'i'-, since 

it is going to be explored by a BABAR Dalitz analysis of 
B- D°K-, D° K°K+K- decays. Excluding the 
channel from the present analysis will allow the re- 
sults of both studies to be more easily combined in the 
future. 

We optimize our event selection to minimize the sta- 
tistical error on the B~-^D°K~ signal yield, determined 
for each decay channel using simulated signal and 
background events. We reject a candidate track if its 
Cherenkov angle does not agree within four standard de- 
viations (cr) with either the pion or kaon hypothesis [14], 
or if it is identified as an electron by the DCH and 
the EMC. Neutral pions are reconstructed by combining 
pairs of photon candidates with energy deposits larger 
than 30 MeV that arc not matched to charged tracks. 
The photon pair invariant mass is required to be in the 



range 115-150 McV/c^ and the total tt*^ energy must be 
greater than 200 MeV in the laboratory frame. To im- 
prove momentum resolution, the invariant mass of the 
two photons from candidate tt'^'s is constrained to the 
nominal tt*^ mass [14]. Neutral kaons are reconstructed 
from pairs of oppositely charged tracks with invariant 
mass within 7.8 McV/c^ 3 cr) of the nominal K'^ mass. 
The ratio between the candidate K'^ flight length and its 
uncertainty must be greater than 2. The w mesons are 
reconstructed from tt+tt^tt'^ combinations with invariant 
mass in the range 0.763 < M(7r+7r-7r°) < 0.799 GeV/c^. 
We define ^jv as the angle between the normal to the co 
decay plane and the momentum in the uj rest frame, 
and Otttt as the angle between the flight direction of one of 
the three pions in the co rest frame and the flight direction 
of one of the other two pions in the two-pion rest frame. 
The quantities cos9n and cos^^^r follow cos^ 9n and 
sin^ Otttt distributions for the signal and are almost flat 
for wrongly reconstructed or false u candidates. We re- 
quire the product cos-^ 9n sin^ > 0.08. The invariant 
mass of a candidate M{D^) must be within 2.5 a of 
the mean fitted mass, with cr ranging from 4 to 20 MeV/c^ 
depending on the D° decay mode. To improve the 
momentum resolution, the candidate invariant mass is 
then constrained to the nominal D'^ mass [14] for all 
decay channels. For D*'^7r^7r+, the invariant mass of 
the (/i~7r+) system, where 7r+ is the pion from the 
and h~ is the prompt track from B~ taken with the kaon 
mass hypothesis [14], must be greater than 1.9 GeV/c^ 
to reject background from B~^D^7r~ , D^^K~n~^ and 
B-^K*°Tr-, K*°^K-n+ decays. We reconstruct B 
meson candidates by combining a D° candidate with a 
track h. For the D'^^K~tt^ mode, the charge of the 
track h must match that of the kaon from the meson 
decay, selecting b ^ c mediated B decays. 

We select B meson candidates using the energy diflfer- 
ence AE = Eg — E*^/2 and the beam-energy-substituted 
mass TTT-ES = V {Eee/2 + pee • pb)'^/E^^ - Pb^ where the 
subscripts ee and B refer to the initial e+e~ system and 
the B candidate, respectively, and the asterisk denotes 
the e+e~ center-of-mass (CM) frame. The ttzes distri- 
butions for B~^D°h~ signals are Gaussian functions 
centered at the B mass with a resolution of 2.6 MeV/c^, 
and do not depend on the D'^ decay mode or on the 
nature of the prompt track. In contrast, the AE dis- 
tributions depend on the mass assigned to the prompt 
track. We evaluate AE with the kaon mass hypothesis 
so that the peaks of the distributions are centered near 
zero for B~^D^K~ events and shifted by approximately 
50 MeV for B~^D°it~ events. The AE resolution de- 
pends on the momentum resolutions of the meson and 
the prompt track and is typically 16 MeV for all D° 
decay modes under study. All B candidates are selected 
with mES within 2.5 cr of the mean value and with AE 
in the range -0.15 < AE < 0.20 GcV. 

To reduce background from e+e^ qq events (with 



q = u,d,s,c), denoted qq in the following, we construct 
a linear Fisher discriminant [15] based on the four event- 
shape quantities L^'^^, \cos9^\,\ cos6%\ and R^'^^. The 
ratio L^^^ between L2 = Pi cos^ 9i and Lq = pi 
is evaluated in the CM frame, where the Pi are the mo- 
menta of charged tracks and neutral clusters not used to 
reconstruct the B (i.e., the rest of the event, ROE), and 
the 9i are their angles with respect to the thrust axis of 
the B candidate's decay products. The angle 6^ is mea- 
sured between the thrust axis of the B candidate's decay 
products and the beam axis, and is evaluated in the CM 
frame. The angle 6*^ is measured between the B candi- 
date momentum and the beam axis, again evaluated in 
the CM frame. The ratio of the Fox- Wolfram mo- 

ments H2 and Hq, is computed using tracks and photons 
in the ROE [16]. The efficiency of the requirement on the 
value of the Fisher discriminant ranges from 74% to 78% 
for B-^D^K-s\gn&\ events and from 17% to 23% for qq 
background events. For the K-k channel, the values are 
87% for signal and 42% for background events. 

For events with multiple S~— >Z)'^/i~candidates (0.4%- 
7.7% of the selected events, depending on the decay 
mode), we choose the B candidate with the smallest = 
Y.^{Mc - {Mc)f/{a'ii^ + Tl) formed from the measured 
and true masses of the composite candidates c, Mc and 
(Mc), scaled by the resolution ctm^ and width Fc of the 
reconstructed mass distributions. Composite candidates 
considered are the B candidate itself (mEs)j t^'^ , and 
UJ candidates. Also F^^ is the only non-negligible width. 

The total reconstruction efficiencies, based on simu- 
lated B D^K events, are 36% {K-tt+), 29% {K-K+), 
29% (7r-7r+), 15% (if°7rO), and 6% {Kluj). 

The main contributions to the background from BB 
events come from the processes B~^D*h~ , B~^D^p~ , 
misreconstructed B~—>D'^h~, and from charmless B 
decays to the same final state as the signal: for in- 
stance, the process B~^K~ K~^K~ is a background 
for B-^D°K-, D'^^K-K+. These charmless back- 
grounds have similar AE and rriEs distributions as the 
D'^K~ signal and are referred to in the following as peak- 
ing BB backgrounds (S" X1X2K-). 

We determine the signal and background yields 
for each decay mode independently from a two- 
dimensional extended unbinned maximum-likelihood fit 
to the selected data events. The fit is performed simulta- 
neously on the B~^ and B~ subsamples. The input vari- 
ables to the fit are AE and the Cherenkov angle 9c of the 
prompt track as measured by the DIRC. The extended 
likelihood C for N candidates is given by the product of 
the probabilities for each individual candidate i and a 
Poisson factor: 

^= j^l[Vi{AE,9c). (6) 

The probability Vi is the sum of the signal and back- 



ground terms, 

v,{\E,ec) = Y.j^'PiE,iK,i ' (7) 
J 

where J denotes the seven signal and background hy- 
potheses D°h, qq{h), BB{h), and X1X2K. N' is the 
total event yield estimated by the fit, and Nj is the 
event yield in each category. We fit directly for the ratios 
R' = R(±) and asymmetries Acp± , as appropriate to the 
decay mode; they enter Eq. (7) through 

Ndo^± = i(lTAg;-)iV^o, , (8) 
Ndok± = ^{IT Acp) Nno^R' , (9) 

where Nijo^^ — N]jo^+ + Ndo^^- and A^p^ is defined 
analogously to Eq. 2. 

The AE distribution for B^->D^K^ signal IS param- 
eterized with a double Gaussian function. The fraction 
of the wide component of the signal shape, its offset 
from the narrow component and the ratio between the 
widths of the two components are fixed to values obtained 
from simulation. The AE probability density function 
(PDF) for B^^D^TT^ is the same as the B^^D"K^ 
one, but with an additional shift, Ai?shift) which arises 
from the wrong mass assignment to the prompt track. 
The shift is computed event by event as a function of 
the prompt track momentum p and a Lorentz factor 
7PEP-II = -E'ee/-£'ee characterizing the boost to the e+e~ 
CM frame: 

A^shift = 7PEP-II [\Jm\+p^ - +p2^ . (10) 

The AE distributions for the continuum background are 
parameterized with a straight line. The AE distribution 
for the BB background is empirically parametrized with 
a Gaussian peak with an exponential tail [17]. The pa- 
rameters of the background shapes are determined from 
simulated events (BB) and off-resonance data (17^) and 
are fixed in the fit. The number of peaking background 
events Nx^x-iK is fixed to values obtained from a study 
of the mass sidebands. The particle identification 
PDF is a double Gaussian as a function of which is 
the difference between the measured Cherenkov angle 9c 
and its expected value for a given mass hypothesis, di- 
vided by the estimated error. The PID shape parameters 
are obtained from simulation. To summarize, the float- 
ing parameters in each of the five the fits are the D^K 
and -D^TT signal yield asymmetries, the total number of 
signal events in D^tt, the appropriate ratios R and R±, 
eight background yields (one for each charge), and two 
parameters of the AE signal shape (common for positive 
and negative samples). 

The results of the fits, expressed in terms of signal 
yields, are summarized in Table I. Figure 1 shows the 



distributions of AE for the ii'~7r+, CP+ and CP— modes 
after enhancing the B~ K~Y)\n:\ty by requiring that 
the prompt track be consistent with the kaon hypothe- 
sis. This requirement is 88% (1%) efficient for h~ = K~ 
{h-=n-). 



TABLE I: Uncorrected yields as obtained from the maximum 
likelihood fit. The quoted uncertainties are statistical. 



D° 


CP 


N{DiT+) 




N{DK+) N{DK-) 


K-TV+ 




12745 ± 120 


12338 ± 120 


954 ±36 918 ±36 


K'K+ 


+ 


1109 ± 36 


1051 ± 35 


51 ±10 113 ±13 


7r~7r"'" 


+ 


390 ± 24 


378 ± 24 


39 ± 9 36 ± 9 






1102 ± 37 


1134 ±38 


100 ± 13 88 ±12 


K°w 




422 ± 24 


403 ± 26 


29 ± 8 18 ± 8 



The ratios R{±), as measured by each fit, are corrected 
to take into account small differences in the selection ef- 
ficiency between B—^DK and B^Dtt. The efficiency 
ratios range from 1.013 ± 0.006 to 1.037 ± 0.010. Their 
uncertainties are due to the statistics of the simulated 
samples and arc considered in the study of systematic 
uncertainties. In the case of D'^ ^s^i lo^tt^ tt'^ , 

the values of Rqp'^ and A^p'^ need to be corrected 
to take into account a possible dilution from a non- 
resonant CP-even background arising from B~^D°h~ , 
D° K°{Tr~ Tr'^)non-uj decays. There is little infor- 
mation on this background. We estimate the corrections 
using a fit to the lu helicity angle in the selected data 
events and find the correction factors to be 1.12 ± 0.14 
for Af-,p_ and 1.00 ±0.01 for Rqp_. The uncertainties in 
the correction factors arc included in the systematic er- 
rors. After applying all corrections, the quantities R±/R 
and Acp± are computed by means of a weighted average 
over the CP+ and CP— modes. The results for the CP- 
even and CP-odd combinations are reported in Table II. 

Systematic uncertainties in Rcp± and Acp± are listed 
in Table III. The uncertainties on the fitted signal yields 
are due to the imperfect knowledge of the AE and PID 
PDFs and of the peaking background yields, and are eval- 
uated in test fits by varying the parameters of the PDFs 
and the peaking background yields by ±la and taking 
the difference in the fit results. A possible ±20% CP 
asymmetry in the peaking background is considered in 
the same way. In the K'^oj channel we also take into ac- 
count the uncertainties in the correction factors due to 
the CP-even backgrounds from /f3(7r^7r+7r'^)non-w 
decays. A possible bias in the measured Acp± comes 
from an intrinsic detector charge asymmetry due to 
asymmetries in acceptance or tracking and particle iden- 
tification efficiencies. An upper limit on this bias is ob- 
tained from the measured asymmetries in the processes 
B-^D"h-, D°^K~ 7r+ and B-^D^p^n-, where CP 
violation is expected to be negligible. From the average 
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FIG. 1; Distributions of AE for events enhanced in 
D°K^ signal: a) B~ D%p+K- ; b) B+ 



c) B- D%p_K-- d) B+ D°cp-K+- B± ^ D"^ , 
D''— >_K'^7r^ with (e) and without (f) signal enhancement. 
Blue (continuous) curve: projection of the full PDF of the 
maximum likelihood fit. Red (long-dashed): B± D°K^ 
signal on all backgrounds. Brown (short-dashed): peaking 
component on qq and BB background. Green (dash-dotted): 
qq and BB background. 



asymmetry, —(1.6 ± 0.6)%, we obtain the limit ±2.2% 
for the bias. For the branching fraction ratios Rcp±, 
an additional source of uncertainty is associated with 
the assumption that Rcp± = R±/R. This assump- 
tion holds only if the magnitude of the ratio r^r between 
the amplitudes of the B~^D^Tr~ and B^^D^ti^ pro- 
cesses is neglected [18]. r^r is expected to be small: 
^-TT ^ f < 0.012, where A w 0.22 [14] is the sine 
of the Cabibbo angle. This introduces a relative uncer- 
tainty ±2r7r cos (5^ cos 7 on Rcp±, where St^ is the relative 
strong phase between the amplitudes A{B~^D^tt~) and 
AiB-^D^TT-). Since |cos(5^cos7| < 1 and < 0.012, 
we assign a relative uncertainty ±2.4% to Rcp±, which 
is completely anti-correlated between Rcp+ and Rcp-- 
We quote the measurements in terms of x± and r^, 

x+ = -0.09±0.05(stat)±0.02(syst), (11) 
a;_ = +0.10±0.05(stat)±0.03(syst), (12) 



= +0.05±0.07(stat)±0.03(syst). 



(13) 



The correlations between the different sources of system- 
atic errors, when non-negligible, are considered when cal- 
culating x± and r^. The measured values of x± are 
consistent with those found from B^^D'^K^ , D'^ 
K'^tt^tt^ decays, and the precision is comparable [11]. 



TABLE II: Measured ratios Rcp± and Acp± for CP-even 
{CP+) and CP-odd (CP—) D decay modes. The first error is 
statistical; the second is systematic. 



D° mode 


Rcp 


Acp 


CP-f 


1.06 ±0.10 ±0.05 


0.27 ±0.09 ±0.04 


CP- 


1.03 ±0.10 ±0.05 


-0.09 ±0.09 ±0.02 



TABLE III: Systematic uncertainties on the observables 
Rcp± and Acp± in absolute terms. 



source 


ARcp+ 


ARcp- 


AAcp+ 


AAcp- 


fixed fit parameters 


0.036 


0.019 


0.010 


0.002 


peaking background 


0.029 


0.037 


0.031 


0.003 


detector charge asym. 






0.022 


0.022 


opp. CP bkg. in K"sLJ 




0.002 




0.007 


Rcp± vs. R± 


0.026 


0.025 






K/tt efficiency 


0.002 


0.007 






total 


0.053 


0.049 


0.039 


0.023 



In conclusion, we have reconstructed B~ -^D^K~ de- 
cays with mesons decaying to non-CP, CP-even and 
CP-odd eigenstates. The combined uncertainties we find 
for Acp± {Rcp±) are smaller by a factor of 0.7 (0.9) 
and 0.6 (0.6) than the previous BaBAR [9] and Belle [10] 



measurements, respectively. We find Acp+ to deviate by 
2.8 standard deviations from zero. We express the re- 
sults in terms of the Cartesian coordinates x± and 
(Eqs. 4, 5). These measurements, combined with the 
existing measurements from B~ —>D'^ K~ decays, will im- 
prove our knowledge of the angle 7 and the parameter 
r. 
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